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J, The discovery of the accelerated expansion of the universe using Type la su- 

"+^ pernovae (SNe la) has stimulated a tremendous amount of interest in the use of 

f—i SNe Type la events as standard cosmological candles, and as a probe of the fun- 

I I damental physics of dark energy. Recent observations of SNe la have indicated 

a significant population difference depending on the host galaxy. These obser- 

T~H vational findings are consistent with SNe la Ni-56 production in star-forming 

^ spiral galaxies some 0.1 solar masses higher — and therefore more luminous — 

^P than in elliptical galaxies. We present recent full-star, 3D simulations of Type 

^N| la supernovae which may help explain the nature of this systematic variation 

C ^ ^ in SNc la luminosities, as well as the nature of the la explosion mechanism. 

C ^ ^ These insights may in turn eventually shed light on the mystery of dark energy 

• itself. 

O Keyword.s: Type la supernovae; computational astrophysics. 



1. Introduction : Observational Properties of Type la 
Supernovae 

The discovery of the accelerated expansion of the universe using Type la 
supernovae^ '^ has stimulated a tremendous amount of interest in the use 
of SNe Type la events as standard cosmological candles, and as a probe of 
the fundamental physics of dark energy. 

Supernovae come in different types, classified by their spectra. They 
are grouped into two main categories based on the presence or absence 
of Balmer absorption lines, indicating the presence or absence of atomic 
hydrogen. Type II contain hydrogen, whereas Type I lack hydrogen. These 



May 4, 2010 0:9 WSPC - Proceedings Trim Size: 9in x 6in FisherlCGAOQ 



types are further subdivided by absorption line features and light curve 
shapes, see Table 1. In particular, Type la exhibit an absorption line of 
singly ionized silicon (Si"*") at 615 nm at peak light. Type Ib/c are further 
classified according to the absence of Si II lines, and the presence or absence 
of He lines, respectively. 

Accordingly, Type la properties are consistent with white dwarf pro- 
genitors. The two leading models include white dwarfs in binary systems - 
either with a non-degenerate companion, the so-called "single-degenerate" 
scenario, or with another white dwarf - the so-called "double-degenerate" 
scenario. While both types of events are likely to occur in nature, and 
contribute to the diversity of explosion energies, there is still considerable 
debate over their relative frequencies, and even what the explosion mecha- 
nism is in each case. We outline some key theoretical and mounting observa- 
tional evidence below based on nearby historic supernovae, which supports 
the identification of typical brightness "Branch normal" supernovae with 
single-degenerate events. 

Table 1. Classification of Supernovae 



Mass — > 




< — Age 


SNIa 


SN II (Iln, IIL, IIP, lib) 


SN Ib/c 



No Hydrogen Has Hydrogen No Hydrogen 

(Si+ absorption) (no Si+) 

White dwarf Progenitor Core collapse of Core collapse 

a massive star (outer layers stripped by winds) 



Tying the variable luminosities of Type la explosions to additional pa- 
rameters is of key significance in any attempt to calibrate the Phillips 
relation. Analysis of light spectra demonstrates that Type la supernovae 
produce a relatively constant combined yield of stable nuclear statistical 
equilibrium nuclei (NSE) of '"'^Ni and ^''Fe and intermediate mass elements 
(IME) Si-Ca of 1.05 ± .09 M©.^ The approximately fixed burned mass 
over a wide range of SNe la luminosities suggests that SNe la share a 
single common explosion mechanism. However, recent observations of SNe 
Type la have indicated a significant population difference depending on 
the host galaxy. In particular, observations have found that SNe la in star- 
forming galaxies decline more slowly than those in ellipticals.* These find- 
ings arc consistent with la SNe '"'^Ni production in star-forming galaxies 
some ^ 0.1 Mq higher than that in ellipticals.^ 

Previous work has focused on the possible influence of the metallicity 
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as a second parameter in explaining the luminosity variance of Type la 
events. ^"^ These results suggest that while the metallicity effect may in- 
deed contribute to the observed variance in Type la luminosities, alone 
it accounts for only a portion of the total variance. ^'^ Consequently, we 
must look elsewhere to other physical effects - including both the central 
density and angular momentum profile of the white dwarf progenitors - 
which may explain the majority of the observed la luminosity variance. 
These progenitor properties may in turn be directly influenced by the sur- 
rounding environment of la event, in particular the accretion rate from the 
companion star. 

To put the single- and double-degenerate white dwarf progenitor mod- 
els into their proper astrophysical context, we briefly consider the obser- 
vational evidence. A number of recent papers have lended strong support 
for the general viability of the single-degenerate channel for the origin of 
Branch normal SNe Type la, in which a progenitor white dwarf accretes 
material from a non-degenerate companion star.^'^" Analysis of the spec- 
trum from the Type la SN 2006X in Virgo suggests that the blast wave 
moved through the circumstellar medium and collided with the ejecta from 
an otherwise undetected red giant companion star.^^ Furthermore, a candi- 
date G-type companion to SN1572 (Tycho) has been identified. ^^ SN1572 
has independently been identified as a Branch-normal la event based on 
analysis of historical light curves and light echo spectra. ^^'^'* 

However, despite this string of recent successes, the single degenerate 
model faces significant theoretical challenges. In particular, models sug- 
gest that a non-rotating white dwarf can burn stably only over a relatively 
narrow mass accretion rate range. ^"^ The addition of differential rotation 
broadens this range. ^^ However, when one introduces rotation to the white 
dwarf structure, one must make a number of assumptions regarding accre- 
tion and internal shear in order to understand the evolution of the angular 
momentum distribution of the white dwarf progenitor. The accretion rate 
in turn sets the central density of the white dwarf at ignition. 

As a consequence of these challenges faced by theoretical descriptions of 
the single degenerate model, and the limited direct observational evidence 
constraining progenitors, there is a significant degree of uncertainty in the 
characterization of the degenerate progenitor. Here we focus on the role 
which the central density (or equivalently, the total mass) of a non-rotating 
white dwarf progenitor plays in determining the luminosity of Type la SNe. 
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2. Physics of Type la Supernovae 

The energetics of the single-degenerate model of Type la supernovae can be 
simply estimated using nothing more than elementary physics. The internal 
energy of N electrons in a fully-degenerate white dwarf is N times the 
characteristic Fermi energy Ep. The Fermi energy can itself be estimated 
for a relativistic electron with momentum p as 

Ef = PC (1) 

Applying the Hcisenberg Uncertainty Principle, p ^ hn^ , where n is the 
number density of electrons. Therefore, 

^ 1 hNTic 

Ep - hn^c (2) 

R 

Where in the second step we have estimated the mean number density 
within a white dwarf of radius R. Including the gravitational binding energy, 
the white dwarf has a total energy of 

E^NEp + Eg^ -^ ^ (3) 

here Eq is the binding energy, and nip is the proton mass. Note that we 
focus here on the essential physics, and have therefore suppressed all factors 
of order unity, and neglected composition effects. 

This elementary analysis reveals a remarkable feature of a fully degener- 
ate star. In a fully degenerate gas, the degeneracy pressure is fundamentally 
independent of temperature. This immediately leads to the result that both 
the internal energy term and the gravitational binding energy term scale 
inversely with radius. Consequently, a fully-degenerate white dwarf cannot 
seek a lower total energy state by an adiabatic spherical compression or ex- 
pansion. This stands in sharp contrast to stars supported by ordinary gas 
pressure, which can indeed minimize their total energy by spherical adia- 
batic compression or expansion. We are led to conclude that the stability 
of the white dwarf is therefore set solely by the sign of the total energy E. 
The critical case is where the total energy E equals zero; solving for the 
maximum total mass Mchandra of the star then yields 

3 2 

-M^Chandra ^ 7=7 I -^ ^ ^Tlpianck ^^^^ I ^1-5 Mq (4) 

\G J ml \ nip J 

This analysis shows the critical mass, which is known as the Chandrasekhar 
limit, to be a combination of fundamental physical constants. Indeed, it is 
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fundamentally set by the Planck mass times a large dimensionless num- 
ber, which is the ratio of the Planck mass to the proton mass, squared. 
This simple analysis, which neglects composition effects, demonstrates the 
Chandrasekhar mass is of order a solar mass; a more precise calculation 
demonstrates it to be 1.4 Mq for a predominantly C/0 white dwarf. 

We take the progenitor to be a carbon-oxygen white dwarf. The mass, 
which is close to the Chandrasekhar limit, undergoes carbon burning, re- 
leasing roughly 10^* ergs/g. The characteristic nuclear energy available for 
a Type la is about 3 • 10^^ ergs ~ 3 foe, where 1 foe is a convenient unit rep- 
resenting 10^^ ergs, a typical Type la luminosity. Therefore, there is more 
than sufficient nuclear energy within the white dwarf progenitor to repre- 
sent not only typical Type la events, but also the more luminous events 
observed. 

In the single-degenerate model of a Type la supernova event, the system 
begins as a binary pair of main sequence stars. As time progresses, the more 
massive star evolves into a giant, accreting gas onto its companion, which 
expands and is eventually engulfed. The core of the giant and its companion 
begin to spiral inward inside a common envelope. Tidal torques cause the 
envelope to be ejected, and the binary separation to decrease. The core of 
the giant then collapses into a white dwarf that begins to accrete gas from 
its aging companion. The white dwarf continues to accrete until reaching a 
critical mass close to but not equal to the Chandrasekhar limit. At a critical 
mass set by the accretion rate from the companion, the central core of the 
white dwarf ignites a nuclear flame that causes the supernova explosion, 
and ejects the companion from the system. 

The precise initial conditions leading to ignition are still poorly under- 
stood. The problem arises because of a large dynamic range between the 
long convection phase (on the order of a hundred years), and the brief defla- 
gration/detonation phase following flame ignition, which is on the order of 
one to several seconds. ^^ At some point during convection, a runaway nu- 
clear burning occurs during which the burning exceeds the neutrino cooling. 
The runaway reaction ignites a deflagration flame slightly off-center from 
the white dwarf, giving rise to a buoyancy force. This buoyant flame bubble 
undergoes subsonic burning while rising toward the surface. The possibility 
of the bubble burning though the entire star in a pure deflagration^^ has 
generally lost support because it tends to burn incfflciently, leaving behind 
a significant amount of nuclear fuel, and producing events which are gen- 
erally underluminous with respect to typical Branch normal la events. It 
is most likely that the bubble cither undergoes a deflagration to detona- 
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tion transition (DDT)^* or leads to a gravitationally-confined detonation 
(GCD).^^ In the case of a GCD, the deflagration bubble breaks through 
surface of the star. Ash is launched into the atmosphere and wraps around 
the surface of the star under the force of gravity. The ash quickly reaches 
the opposite side of the star and collides with itself. This collision creates 
a supersonic detonation front that unbinds the star. 

At present, both the DDT and GCD mechanisms have their advan- 
tages and disadvantages. The DDT model can be tuned to be consistent 
with observations. However, in the absence of a full first-principles under- 
standing of the detonation transition, DDT simulations set the detonation 
transition as a parameter. The current belief is that detonation can occur 
when the flame is ripped apart by turbulence, and transitions into the dis- 
tributed burning regime. However, the precise conditions under which this 
transition occurs remain a matter of intense research, so that the transition 
density is largely a free parameter in the simulations. On the other hand, 
the GCD mechanism can successfully produce detonations without fine- 
tuning, though initial models typically overproduced ^®Ni, underproduced 
intermediate mass elements (IME), and generated overluminuous events. 
We address one aspect of the luminosity problem in the following section, 
by varying the central density of the white dwarf progenitor models in the 
simulations. 

3. Simulations of Type la Supernovae 

We begin with a carbon-oxygen white dwarf with a pre- ignited flame bubble 
of initialized size and location. The model is run through detonation. The 
evolution of the white dwarf and bubble has a broad range of length scales. 
On the largest scales, one must be able to follow the expansion of the 
supernovae out to several tens of thousands of kilometers to capture the 
homologous expansion phase, and on the smallest scales, one must begin 
to capture the flame physics, which extends down to the laminar flame 
thickness on centimeter scales. We employ the Paramesh library within the 
FLASH code to implement an adaptive mesh refinement (AMR) mesh. Even 
the power of AMR still does not allow one to avoid the enormous dynamic 
range between the large-scale physics of the explosion and the flame scale 
- over 10^ in linear dimension - so we incorporate a thickened model of the 
flame surface, which artificially thickens the fiame over ~ 4 grid cells. 

We evolve the simulation using a comprehensive multiphysics model, in- 
cluding the coupled equations of hydrodynamics, self-gravity, and nuclear 
combustion. Specifically, we incorporate the Euler equations of hydrody- 
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namics : 



dt 

dpv 

1h' 

dpE 



V • (vp) ^ (5) 

- V- (w«)pw) = -VP-pV$ (6) 

+ V • [v{pE + P)] = pw • V$ + penuc (7) 



with source terms for self-gravity and nuclear energy release. Here p is mass 
density, v is velocity, P is pressure, $ is gravitational potential, e„„c is the 
specific nuclear energy release. E is the total energy density, given by 

E = p(^U+^v'^ (8) 

Where U is the specific internal energy. 

The Euler equations are coupled to Poisson's equation for self-gravity 

V2$ = AnGp (9) 

and an advection-difFusion reaction model of the thickened combustion front 

^+^.y^ = kV^^ + -R((j)) (10) 

Ot T 

R{cj,)^^{<f>-eo){l-cj, + e,) (11) 

Here (f> is a scalar progress variable which monitors the advancement of 
the flame surface, and sets the nuclear energy release function Cnuc above. 
K is a parameter representing the diffusivity of the thickened flame, r is 
the reaction timescale, R is the reaction term, which is set by a lowered 
Kolmogorov-Petrovski-Piscounov (KPP) binomial, specifled by the three 
parameters, /, eo, and ei. 

Our first sucessful 3-D simulation of a Type la detonation was a cold 
white dwarf model in initial equilibrium with initial mass 1.36 Mq. The 
nuclear bubble was ignited within a spherical region slightly offset from 
the center of the white dwarf. The supercomputer simulation "marches" 
this condition forward in time in 3-D, using full equations describing the 
flame, hydrodynamics, and self-gravity. An numerical analysis yields a crit- 
ical conditions for initiation.^" We found that the critical conditions for 
degenerate white dwarf matter are robustly met in our 3-D simulations. 
We have confirmed that detonation arises independent of the resolution in 
the detonation region, and also for a wide variety of initial bubble sizes and 
offsets. Simulations of the GCD model produce intermediate mass elements 
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Fig. 1. Density versus radius for a series of cold white dwarf models, ranging from 
1.355 Mq to 1.385 Mq. The plot focuses on the innermost region of the white dwarf to 
emphasize the sensitive dependence of central density on total mass. 



at a velocity coordinate ^11, 000 km/s, creating a layered structure of IME 
and Fe peak (NSE) products similar to observation.'^ 

However, despite these successes, initial simulated 3D GCD models gen- 
erally underproduced intermediate mass elements and overproduced Ni. 
Consequently, the models are generally too luminous in comparison to 
Branch normal la events. The underabundance of IME suggest that ad- 
ditional pre-expansion is required to produce Ni and IME abundances con- 
sistent with observation. Moreover, other 3D simulations were found to 
pre-expand significantly and lead to failed explosions. ^^ 

One factor influencing this outcome is the uncertainty in the central 
density of the progenitor model, which is fundamentally set by the accre- 
tion rate in the single degenerate scenario. As shown in figure 1, a small 
change in the total mass of a cold 50/50 C/0 white dwarf progenitor model 
from 1.355 Mq to 1.385 Mq leads to more than a factor of four change in 
central density : from 10^ gm/cm'^ to 4.4 • 10^ gm/cm^. This increase in 
central density in turn directly impacts the nuclear energy release during 
the deflagration phase. Fundamentally, this is due to two reasons. First, the 
critical wavelength for Rayleigh- Taylor instability Ac ex: Sf /g is sensitively 
dependent on the density through the dependence on the laminar speed 
5*;. Here g is the gravitational acceleration. Higher central densities lead 
to a larger critical wavelength for the development of the Rayleigh- Taylor 
instability, which in turn suppresses the growth of turbulence during the 
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early stages of the burning, and leads to an enhancement in the burnt mass. 
Second, the increased total mass of the progenitor with increased central 
density causes it to be closer to the Chandrasekhar mass, and therefore 
more easily pre-expanded than models with lower central density. 

The net effect on the energy release during the deflagration phase in the 
GCD model is dramatic, as shown in figure 2. An increase in total mass 
by just one percent, and an increase in central density by a factor of two, 
leads to a four-fold increase in the total fractional energy release during the 
deflagration phase. This leads to a significantly enhanced pre-expansion 
rate in the case of the higher central density model. 
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Fig. 2. A plot of the nuclear energy released as a function of time, shown for a range 
of stellar masses, from 1.365 Mq to 1.385 Mq. 



4. Conclusion 

The GCD model simulations demonstrate that successful self-consistent 
detonations of Type la supernovae in 3-D are possible without having to 
be initiated by hand. The nuclear energy release during the deflagration 
phase is found to be sensitively dependent upon the central density of the 
white dwarf progenitor model. With a greater pre-expansion afforded by the 
greater nuclear energy release in higher central density models, the GCD 
mechanism can lead to observed Branch normal luminosities, and typical 
levels of intermediate mass elements (Si -Ca). 
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